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ABSTRACT 


The design of a supersonic missile seeker lens system 
uSing an aerodynamically efficient lens, a scanning mirror 
and detector was accomplished through the use of ray- 
tracing routines. The design of the first element lens 
uses a gradient in refractive index to overcome the severe 
handicap that the pointed shape imposes. This proolem has 
been previously solved for positive orientations of the 
center of symmetry of the Gradient Refractive Index (GRIN) 
material with respect to the lens. In this study the 
negative orientations are resolved and the lens studies for 
a range of orientations and strengths of gradient. The 
best lens is then used in conjunction with the mirror and 
detector in a ray-tracing scheme. Such a system was found 


to be feasible and worthy of further study. fh —_ 
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I. INTRODUCTION 


A. OUTLINE OF LENS DESIGN TO DATE 

The design of a 'pointed' aerodynamically efficient 
missile seeker lens has been the topic of a number of 
studies. Frazier [Ref. 1], and Terrell [Ref. 2], first 
studied the idea of using a pointed lens as a means of 
obtaining some degree of aerodynamic efficiency. They also 
introduced the idea of using Gradient Refractive Index 
(GRIN) material in order to overcome the difficulties in 
imaging imposed by the odd shape. Amichai [Ref. 3], and 
Carr [Ref. 4], studied the GRIN concept in more detail, 
Amichai by defining the outer surface as a cone and solving 
for the inner surface and Carr by defining the inner 
surface as a cone and solving for the outer. The main 
reference for this author's study was drawn from Carr's 
work and a review of that thesis would be helpful if the 
reader seeks a more detailed background from that given in 


the next few pages. 


B. GRADIENT REFRACTIVE INDEX (GRIN) 

Although the theory for GRIN has been developed and 
Known for some time it was not until fairly recently that 
materials have become available that allow the theory to be 
put into practice. A gradient in index will cause light 
rays to bend in an otherwise homogeneous material. Thus, 
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the lens designer now has a new dimension of flexibility 
which will allow him to solve problems that could not be 

4 formerly solved. A proper selection of the manner in which 

the gradient varies is also important. The designer would 

normally select certain degrees of symmetry. Three are 

_| most common; axial, radial or spherical gradients. 

GRIN material in the shape of rods now commonly use 
axial and radial gradients in such devices as fiber-optic 
cables and photocopying machines. Spherical gradients are 
not available because spherical gradient material has not 
yet been produced. Spherical gradients do offer an 


advantage in that they lend themselves more easily to 


a ee 


mathematical description since closed form solutions for 


the ray paths can be obtained. 


a A good introduction to gradients and the materials 

| currently available for production of large scale GRIN lens 
systems is provided by Moore [Ref. 6], and Light [Ref. 7]. 
These also give the reader an idea of the advantages that 
may be realized along with the problems associated witn the 


design and production of such items. 


C. AIM AND INTENT OF THIS STUDY 

The aim of this study is to complete the design of the 
e Gradient Index Seeker Lens (GISL) initiated by Carr. This 
design uses a spherical gradient and 'pointed' lens shape. 


This study does not determine the aerodynamic efficiency of 
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the design but rather ensures that the lens has a 
generally pointed character that will also allow some 
degree of imaging for off-axis objects. The spherical 
gradient is chosen since it Seems to fit the geometrical 
shape of the lens and also since the three-dimensional 
analysis of rays is much simpler. 

This design is accomplished through the use of ray 
tracing techniques. Thus, no actual study of the lens 
material nor the manner in which such a lens could be 
manufactured is intended. It is the hope of the author 
that studies of this nature will provide the impetus for 
the development of suitable GRIN material with spherical 
gradients. 

Once the lens design was completed a 'best' lens was 
chosen and a simple detector system applied to it to see if 


precise definition of off-axis targets is possible. 
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II. GRIN LENS DESIGN 


A. GRIN LENS DESIGN METHODOLOGY 

The method by which the GRIN Lens was designed in 
Carr's thesis [Ref. 4], is the subject of this section. A 
cross-section view of the top half of the lens is shown in 
Figure 2.1. 

As can be seen, the lens is designed by defining the 
inner surface as a cone and then constructing the outer 
surface so that a backwards ray trace froin a point source 
will result in parallel lines emerging. This is done by 
tracing rays from the back edge of the lens towards the 
front surface and finding the point of intersection of the 
cay and front surface line. The front surface line is 
found by using the previous ray solution as described 
below. The intersection of the ray and the line is 
determined using an iteration procedure. Once a solution 
is found the front surface line that will cause the ray to 
emerge parallel with all other rays is determined. This 
front surface line is then used to find a point solution 
for the next ray. avout 1000 rays are traced starting at 
the outside edge of the inner cone and working inwards 


towards the center. Details of the iteration process is 


included in later sections. 


Once the lens outside surface has been constructed the 
lens is then studied to see how well it will image three 
dimensional skew rays that enter the lens at varying 
degrees of incidence. To accomplish this, a grid is placed 
in front of the lens and then rays are traced from the 
intersection points of the grid through the lens and onto 
the image plane. The grid is then tilted at successive 
angles and the resulting 'spots' on the image plane 
recorded for each case. The manner in which these spots 
move On the image plane and their size are the major 


performance characteristics of the lens. 


B. GRIN RAY PATH DESCRIPTION 
Following Marchand [Ref. 5], GRIN ray paths in 
spherical gradients are found through the solution of the 


appropriate line integral: 


a eee rey oe aro U4 (2.1) 
where 
e = nrsin = ngorgsin g (2.2) 
and 
2 4} 
It should be noted that the value of is -1.9 if the 


radius is decreasing with angle and +1.0 if the radius is 


increasing with angle. 
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The choice of index gradient is instrumental in 
allowi:ig the integral to be solved in closed form. In this 
study the index is defined as: 

n2 = (a+ br2/R,2) 
The constant a defines the strength of the index at the 
origin and the constant b defines the strength of the 
gradient and whether it is negative or positive. 

As can be seen in Figure 2.2, the ray path is described 
by the constant e in conjunction with the parameter °, 
the angle between the radial arm of the coordinate system 
and the direction of the ray at any point. 

The solution of the integral in equation 2.1 results in 
the following expression: 


1 


[sin™* (arg) ) - sin “(arg)! (2.3) 


where 
a- 2e/t,° 
arg, = Pe eae Py ae 
0 tat + 4be*/R,*)* Y 


and 


2 
a- 2e/r 
arg = op ee 
as abe 7,6 


(a 


or, expressed in terms of the radius, 


Y2\e| 
{at (a*+4be*/R,*) sin [-2¢e(*3 


. = 1 
9) 7 sin (arg, )]- 


C. GRIN LENS GEOMETRY 

A factor which must be considered is the manner in 
which the lens is cut out of the GRIN material. Figure 2.2 
Shows an orientation of the material with respect to the 
lens such that the center of symmetry of the material is 
inside the lens. Carr has defined this as a negative value 
of OB, the distance from the point O to the point B. 
Carr was able to work out the lens problem for positive 
values of this parameter only, thus the major portion of 
this study is to develop methods that enable solutions for 
negative values of OB. In so doing, two factors have to 
be overcome: 

a) ray matching - When the ray path radius initially 
decreases and then increases with change in angle the ray 
path equation must be matched at the point of minimum 
radius. 

b) near-radial lines - Rays that show very little 
curvature are difficult to describe in terms of radius and 
angle and these cause problems in the iterative solution 


scheme. 
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The solution to the ray matching problem is relatively 
simple and is discussed in the next section. The 
resolution of near-radial lines is a more complex task and 


is covered in the following chapter. 


D. RAY MATCHING 

As stated previously the ray path equation changes 
Signs as the radius of the path decreases and then 
increases. Such a path is shown in Figure 2.3. Simply 
changing sian in the equation as the point of minimun 
radius is passed does not produce the desired result, but 
rather has the eftect shown in Figures 2.3 and 2.4. 

The resolution of this problem is fairly simple. At 
the point of minimum radius the arcsine function in 
equation 2.3 containing this radius (arg) reduces to - /2. 
This is due to the fact that at this point the value of 


is -/2. Thus: 


sin ;=l1 
and 
e2 = nér2 
Thus: 
2 a-2(atbr’/R 2) 
arg = a~2n = PS ae rE ee aL RE ah = 1 
Few athe ay eee ee) 7 ? Ps 
(a*+4bn747/R,*)? 2 (a*+4b(atbr?/R*)x°/R,*) 
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and 
sin! (arg) =-1/2 
If the radius is initially decreasing as Shown in 
Figure 2.3, we can find the angle of the minimum radius. 


We call this the critical angle or 3: 
. 1 ‘ 
{sin (argy) + 1/2! 


The positive and the negative forms of the equation 
should both contain the point of minimum radius, in other 
words, by approaching the point from the left or the right 
the same result should be obtained. The fact that this 
does not happen is shown in Figure 2.3 and reveals the need 
to match the rays at the point of minimum radius. 

This matching is achieved by simply adding a constant 
to either the positive or the negative form of equation 2.3 
and then matching at the point of minimum radius in the 
following manner: 

PR = S8,¢+ Cc 


thus 


[n/2+sin™* (arg, )! = + {x/2+sin” *(arg,)! +c 


= 
2 
The constant c can then be found: 


arene e'n/2tsin + (arg, )! 


This constant must now be applied for all points past 


the point of minimum radius; in our example all points in 
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the increasiny radius portion of the ray must be found 


using the form: 


% 5 'sin”* (arg, ) = sin’ t(arg)} + ¢ 

In the iteration scheme the ray path equation must be 
more carefully handled. First, the critical angle and the 
constant c must be found. Second, the angle of the 


ray must be checked to see if it has passed the critical 


angle. If it has then the constant c is added to 


. equation 2.3 and the point recalculated. The details of 
the iteration methods will be described in the following 
° chapter. 
For further clarification Figure 2.5 shows the regions 
in which ray matching is required for a certain position of 
Cc the center of symmetry. Note also that were the value of ° 
OB positive then all rays in the lens would never reach a : 
minimum radius point. Thus ray matching is a problem 3 
old unique to the negative OB case. a?“ 
‘eo 
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Figure 2.5 Ray Matching Regions. 
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III. NEAR-RADIAL LINES 


A. DESCRIPTION OF DIFFICULTY 

As mentioned previously, the ray paths that are nearly 
radial are difficult to describe mathematically in terms of 
radius and angle. This is due to the simple reason that a 
Straight line originating from the center of the coordinateé 
system has no change in angle and an undefinable chanye in 
radius. 

The region in which this problem occurs in the lens 
design is shown in Figure 3.1. As can be seen, near-radial 
lines mean that the angle ww and thus the constant e 
will be. very small, so small as to cause equation 2.4 to 
approach a singularity and thus break down. In the 
iterative scheme, as the front surface solutions are found 
from the outside edge of the lens inwards, the values of 
and e become increasingly smaller if the center of 
symmetry is inside the lens (OB negative). when this 
region is reached, equation 2.4 reduces the ry _ for all 
angles and futher construction of the front surface is 
rendered impossible. This effect is shown in Figure 3.1. 

Of note is the fact that if the radial line region is 
passed successfully, the rays hegin to curve in the 


Opposite sense, that is they are always convex to the 
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radial lines in the case of a negative gradient and concave 
to the radial lines if a positive gradient is used. 

The resolution of the near-radial lines proved to be a 
difficult task and became a major portion of this study. 
The manner in which it was overcome is the subject of the 


following sections. 


B. PRELIMINARY WORK - IMPROVING THE ITERATION SCHEME 

The first step in the resolution of the problem was to 
Simplify the iteration scheme. In his work Carr developed 
a scheme based on the geometry of the rays and coordinate 
system [Ref. 4], which worked well for all positive values 
of OB (that is, for all orientations of the center of 
symmetry to the left of the inside surface cone. See 
Figure 2.2). It was decided instead to use the Newton 
Raphson approach, which is simpler in concept and not 
dependent on geometry. 

Carr used the Newton-Raphson method in the second 
portion of the design involving the skew ray trace, but he 
did not use it for the lens construction, 

The manner in which this scheme works is now described 
so that the reader will be able to better understand tne 
concept of how all iterative methods work for this design, 
as well as how the ray matching previously described fits 


in. 


31 


Y AXIS 


1. 


1. 


0. 


-5 


0 
-1.0 


Figure 3.1 


OSYMB = -0.5000 
A = 2.2500000 


B = -0.9843750 


front surface showing | 


affect of singularity | 


region of 
near-radia! 


-0.5 -0.0 0.5 


Ray Trace Showing Region of Radial Lines anc 


Effect of Singularity. 


32 


lines 


wy, 
' 


The geometry associated with this plan is shown in 
Figure 3.2 and a flowchart outlining the iteration logic is 
shown as Figure 3.3. As can be seen with reference to the 
two figures, an initial point is first chosen by drawing a 
straight line from the inside surface intercept to the 
Slope of the front surface obtained from a previous ray 
path solution. Newton Raphson seeks to find the intercept 
of this front surface line with the curved ray path. This 
is done by uSing the initial angle to find the radius and 
derivative with respect to angle for each line and then 


calculating a new angle according to the relation: 


ah ee eh fs gr 
“ts ~ “gr 


This new angle is then used to recalculate the radii 
and derivatives and thus an iterative loop if formed. Once 
the difference in angle is small enough the iteration is 
stopped and the solution considered to be found. 

Ray matching is also employed in this scheme. As each 
new angle is calculated, it is checked to see if the 
critical angle has been passed. If it has, the constant ¢ 
of equation 2.5 is added to the angle and thus the proper 
radius and derivative will be then found. 

Once a solution has been found it is used to find a new 


front surface slope by determining that slope which will 
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cause the refracted ray to emerge parallel to the axis. It 


is this slope that is used to find the solution for the 
next ray in the scheme. 

When this new iterative procedure was implemented it 
was found to work efficiently for both the positive and 
negative gradients, but it did not resolve the problem of 
near-radial lines. 

The reader should refer to Appendix A for a listing of 
program GISL, which is the ray tracing program originally 


written by Carr and modified by the author. 


C. PRELIMINARY WORK - DOUBLE PRECISION 

The next avenue that was explored was to place the ray 
tracing scheme in double precision. This was considered 
necessary Since very small differences in angle were 
required in the region of the radial lines. 

Once this was installed in the routine it was found to 
improve significantly the accuracy of the calculations and 
succeeded in delaying the point at which the iterations 
broke down. It did not solve the crux of the difficulty 


however. 


D. PRELIMINARY WORK —- CHECKING ACCURACY 
At this point it was decided to develop a scheme to 
check the accuracy of the front surface solution. The 


manner in which this was done is shown in Figure 3.4 
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The ray trace was completed in a forwarj direction from 
the front surface to the focal point in a two-dimensional 
fashion. These rays were caused to originate from the mid 
point of each set of points that define the front surface. 
The slope of the front surface for each of these rays was 
taxen to be the line drawn between the two adjacent boints 
on the front surface. In this fashion the worst case was 
chosen since the skew ray trace used to determine the 
performance of the lens uses incident rays that can hit the 
front surface at any point. 

These rays were then traced tnrough tne lens and back 
towards the focal point. The x-coordinate of the focal 
point was defined and the y-coordinate calculated from the 
ray trace to see how close it would come to the 
y- coordinate SE the focal point originally used to start 
construction of the front surface. A listing of this 
routine is shown as Nemec B. 

Another benefit of this two-dimensional ray trace was 
that it allowed the author to confirm the methods by which 
a forwards ray trace is done. In this manner ray matching 
and sign convention in angles was checked. Also this ray 
trace is different in that a solution between the ray path 
in the lens and the back surface, rather than the front 
surface, was now required in the iterative solution scheme 
and this could now be verified to ensure that it was 


properly done. 
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Another method of checking the accuracy of the solution 
would be to study the front surface smoothness in tne 
region of near-radial lines. Although not as precise as 
the above method this idea allowed the author to confirm 
that the front surface was sufficiently smooth. By this 
means it was also possible to see a physical representation 


of the lens in the region of the singularity. 


E. FIRST PROPOSED SOLUTION - STRAIGHT LINES 

The first proposed method of solving the mathematical 
difficulty was to simply replace the near-radial lines with 
completely straight lines just before the singularity was 
reached. This could be done by following the decrease of 
the parameter te until it became small enough to ensure 
that the lines were nearly straight. 

This idea was implemented but it was found not to be 
sufficiently accurate. It must be emphasized that an 
inaccuracy in even a single point could not be tolerated 
Since this would mean that all succeeding points in the 
iterative solution scheme would carry that error as well. 

For these reasons the straight line solution was 
rejected and other means of solving the difficulty were 


sought. 


F. SECOND PROPOSED SOLUTION - REDUCED EQUATIONS 
The second method that was used to resolve the problem 
of near-radial lines was to reduce equation 2.4 by 


38 


3 oe 
Sta me ae > eee Gee NER aera 


ee 


start of trace 


end of trace 


Figure 3.4 Geometry of Accuracy Check - Forwards Ray Tracc. 


39 


cme at Rat an a EE a a ee aa ee i ee 


employing small angle aproximations, binomial expansions 
and curve fitting. The details of how this was done 
4 follows. 


Recalling equation 2.3: 


i ae nt 


A = S$. = 


. ml] .-i] 
{ - g)! 
9 sin (arg ) Sin ~ (ary) 


2 
The arguments of these arcsine functions were reworked 
as follows: 


olf a-2(atb") yy? 


arg, = oe er ee 
0 (a*+4b'(atb") 9%)? 2 


: and 
e a-2(atb') 1 7ey2 70? 
arg = —————__———— 
(a*+4b' (atb') yy 2)1/? 
where a 
My co sini, 


b' = bro2/Rz2 
Note the angle +g is very small. 7 
These arguments were then reduced using a binomial 


expansion: 
2 


a(atb")* 1, 


2 
a 


arg, =l- 
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Figure 3.5 Expanded View of Front Surface -- Straight Line 


Method. 


and 


y 2 7 
tyy, Pie. ’ 
4 2(atb dy ( 5) a+b!) 
arg = 1 - chat eee ar ae rare es 
2 
a 
a Again the ty term ensures that the expansion is valid. 


A function of the type arcsin(x) where x is known to 
be close to the value 1 can be curve fitted into the form: 


vi sin-l(x) = 1/2 = [2(1 - x)]1/2 
This eyuation applies to both arcsine functions as the "9 


term ensures that the argument is close to l. 


. 
When placed in the equation the arguments result in the 
following final expression: 
. be sade _—s 
6 = 9.42 [(atb') - [(atb') (2 ato'y]}/7 1 
0 2 re? 
In terms of radius: 
ot 7 (a(atb'))'/? 
BT ee 2 72 
([—- (9 - 5.) = (atb')]* = b'(atb')! 
EY 0 


The reduced equation was then implemented in the 
iterative front surface solution scheme by constantly 


checking it against the complete equation. When the 
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difference between the two became small, the reduced 


equation was subsequently used in solving for the front 
surface. The reduced equation was then turned off when the 
rays had progressed an equal distance on the other side of 
the singularity. See Figure 3.2 for clarification. 

The equation worked well in resolving the front surface 
for the lens construction and did produce a smooth surface 
Similar to that shown in Figure 3.5. 

It did not do well when tne accuracy program check was 
applied. One reason is that although it was easy to 
implement the equation for the lens construction it was not 
so easy to complete the forwards trace since these rays hit 
the front surface randomly making it difficult to know when 
to use the reduced equation. In practice the reduced 
equation could be implemented when the absolute value of 


the parameter became smaller than that used in the 


¥9 
lens construction phase during the time that the reduced 
eguation was used. There was then no way of telling how 
accurate the front surface was, in other words, there was 
no method to show how the rays would physically react in 
the region of the singularity. Thus the reduced equation 


method would work but it is not an entirely satisfactory 


solution. 
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Program CHECK was tried without using the reduced 


equation but disastrous results were obtained for those 
points in the region of the singularity. 
For these reasons it was decided that another means of 


resolving the problem was required. 


G. THIRD PROPOSED SOLUTION - MARCHAND METHOD 

The third proposed solution is drawn from Marchand 
{Ref. 5] p. 27, in which he proposes a means of improving 
the ray trace formula so that nearly radial lines will be 
more easily described. 


The method is based on redefining variables as 


follows: 
9 = eM (3.1) 
where 
Me, ee ess, 
0 Yo r(n*re - e4)!/2 


Marchand describes the ray path using the parameters x 
and 8 such that: 


& = sin?/siniy 
(3.2) 


a = coss - PCOS ig 


This enables the coordinates of the ray to be described 


as: 
x Pp 
x= la + 75) 
0 0 
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Upon substitution of equations 2.2 and 3.2, equation 

3.1 becomes: 
= -NgrtgMg sinc - 

The parameter can now be more accurately solved for 
small angles of uSing a Maclaurin series for the sinc 
function. 

This method was then implemented in program GISL using 
an iteration scheme based on Cartesian rather than polar 
coordinates. 

This method was not found to improve the accuracy of the 
ray trace in the region of near radial lines and was not 
adopted. 

One reason for this may lie in the fact that in our 
solution to the integral in equation 2.1 the value e no 
longer appears at the front of the expression. Thus there 
is no advantage to redefining the variable as shown in 


equation 3.1. 


H. FOURTH PROPOSED SOLUTION — ANGLE IN TERMS OF RADIUS 
The fourth method and the one that proved successful 


involved the solution of the front surface slope and the 
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GRIN ray using the form of the ray trace equation shown in 
equation 2.3 rather than that of equation 2.4. 

The reason for this choice is shown in Figure 3.7. 
This figure shows that when the rays are nearly radial a 
small change in angle will result in a huge change in 
radius. Thus the form of the ray tace equation using r( ) 
is inherently unstable. Previously we have used this form 
in order to sake use of the NewtonRaphson technique. It 
was thought that if the iteration scheme could be reworked 
so that the form i(r) of the equation could be used then 
better results might be obtained. 

A suitable iteration method was found and is 
illustrated in Figures 3.8 and 3.9. 

The figures show that an initial point was chosen in 
the same manner as before, that is by drawing a straight 
line from the ray refracted at the inner surface out to the 
front surface slope. This point provides the radius for 
equation 2.4 so that a new angle can be calculated. This 
angle is then used to find the corresponding radius out to 
the front surface slope. This new radius is then used in 
equation 2.4 to find the next angle and thus an iterative 
loop is formed. 

Thus a simple scheme was found. It should be noted 
however that this method will only converge when the 
initial angle is greater than 90 degrees. The reason for 
this lies in the geometry and is shown in Figure 3.10. 
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Thus the iteration will only converge when the rays are in 
the region of the nearly radial lines. This does not cause 
a difficulty since the Newton Raphson method, modified to 
include ray matching, can find solutions efficiently in the 
region where the initial angle is less than 90 degrees. 

Figure 3.10 is also drawn to show the reader an example 
of a forward ray trace. 

This simple scheme was implemented and found to work 
successfully for both the lens construction and the skew 


ray trace; see [Ref. 8]. 
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Figure 3.9 Angle in Terms of Radius Iteration Geometry. 
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Figure 3.10 Diverging Iteration - Forward Ray Trace. 
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IV. GRADIENT INDEX SEEKER LENS RAY TRACE RESULTS 


A. FRONT SURFACE ACCURACY 

As mentioned previously, the front surface of the 
Gradient Index Seeker Lens was checxed for accuracy by 
completing a ray trace from the midpoint of each set of 
solutions. A sample of the output from program Check is 
shown in Table I so that the reader will have an indication 
of the accuracy of the ray trace in the region of near- 
radial lines. Note that the most radial lines occur when 


the absolute value of the angle is smallest. 


B. SPOT SIZE RESULTS 

The results of program GISL showing the variation of 
spot size are shown in Figure 4.1. This figure shows how 
the spot size varies as gradient strengths and orientations 
of the center of symmetry are changed for a median angle of 
incidence of 0.3 radians. 

Note that all orientations of the center of symmetry 
are now inclvded. The studious reader should refer to 


Carr's work [Ref. 4], to compare results. 


C. SELECTION OF A BEST LENS 
The best lens was selected by determining those 
parameters that would result in a small spot size and also 


“ause the least number of rays to be reflected or refracted 
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Outside the missile diameter. Figuce 4.1 shows that these 


criterion were not easily discernable and a judgemental 
decision was required. 

The reader should note that a value of 2.25 was chosen 
for the parameter a. This was done since the resulting 
value of 1.5 for the index of refraction at the center of 
symmetry is common to a large number of possible lens 
materials. Also of note is the fact that the focal point 
was chosen to be 2.0 units along the x-axis from the inner 
surface cone. This was done since this value resulted in 
smaller spot sizes generally. 

With these criterion in mind a value of -0.2 for 02 
and a 5% negative gradient were chosen as parameters for 
the best lens design. 

A more detailed study of the lens parameters could be 
accomplished with use of an optimization program in which a 
large number of parameters could be varied in order to find 
the best lens, This would achieve a more satisfactory 


result than the one shown. 
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TABLE I 


FRONT SURFACE ACCURACY 


4 | 
{Line Number [ Angle “yg [7 Error in y- number of | . 
| J} (Radians) | coordinate | iterations | . 
| ] a ples vec Sy Sod ei 
| 895 | 0.0411 | 0.0008 | 4 | 
| 896 | 0.0399 | Q .0007 | 4 | 4 
| 8397 | 0.0369 | 0.0008 | 4 | : 
| 898 | 0.0349 | 0.0009 | 4 | 
| 899 | 0.0328 | 0.0007 | 4 | 
| 900 | 0.0308 | 0.0009 | 4 | 

- | 901 | 0.0288 | 0.0008 | 4 | 
| 902 | 0.0268 | 0.0008 | 4 | 
| 903 | 0.0247 | 0.0009 | 4 | 
| 904 | 0.0227 | 0.0009 | 3 | 
| 905 | 0.0207 | 0.0008 | 3 | 
| 906 | 0.0137 | 0.0009 | 3 | 
| 907 | 0.0167 | 0.0008 | 3 | 
| 908 0.0147 | 0.0008 | 3 | 
| 909 | 0.0127 | 0.0010 | 3 | 
| 910 | 0.0107 | 0.0009 | 3 | 
| 911 | 0.0087 | 0.0008 | 3 | s 
| 912 | 0.0067 | 0.0010 | 3 | % 
| 913 | 0.0048 | 0.0009 | 3 | fae 
914 0.0028 | 0.0010 | 2 | 
| 915 | 0.0008 | 0.0010 | 2 | 
| 916 } -0.0010 | 0.0010 | 2 | 
| 917 | -0.0030 | 0.0010 | 3 | 
| 918 | -0.0050 | 0.0010 | 3 | 
919 { -0.0069 | 0.0009 | 3 | 
| 920 | -0.0089 | 0.0011 | 3 | 
| 921 | -0.0108 | 0.0910 | 3 | 
| 922 | -0.0127 | 0.0010 | 3 | 
| 923 | -9.0147 | 0.0009 | 3 | 
| 924 | -0.0166 | 0.0011 | 3 | 
| 925 { -0.0185 | 0.0010 | 3 | 
| 926 | -0.0204 | 0.0911 | 4 | 
| 927 | -0.0223 | 0.0011 | 4 | 
| 928 | -0.0243 | 0.0010 | 4 | 
| 929 | -9.0261 | 0.0012 | 4 | 
| 930 | -0.0281 | 0.0911 | 4 | 
| 931 {| =-0.0300 | 0.0012 | 4 | 
| 932 | -0.0318 | 0.0012 | 4 | 
| 933 | -0.0337 | 0.0012 | 4 | 
| 934 | -0.0356 | 0.0011 | 4 | 
| 935 | -0.0375 | 0.0013 | 4 | 
| | | | | 
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V. MIRROR AND DETECTOR SYSTEM 


A. DESCRIPTION OF SYS‘TEM 

A scanning mirror and a detector were added to “the GRIN 
best lens in order to determine if precise definition of 
off-axis targets would be possible. The design of these 
two new elements is shown in Figure 5.1. 

The mirror is made to ‘nod' at a high rate and thus 
scans the spot across the detector. The detector consists 
of two sections, as shown in Figure 5.1. Thus, when the 
spot is scanned across the two plates a signal from each 
detector will be generated and will result in the form 
Shown in Figure 5.2. This arrangement will allow an 
accurate measurement of the mirror angle for which the spot 
crosses the axis. The electronics of the detector signal 
processor could then be arranged to relate the mirror angle 
to the angle of incident radiation so that the precise 
line-of-sight angle of an off-axis taryet will be known. 

Note that this is a simple detector system that will 
only resolve angles in one plane. A more sopnisticated 
system along the same lines could be designed with a four- 
element detector and a more intricate method of scanning 
the spot across the detector. It was decided to keep the 


mirror and detector system as simple as possible since the 
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aim of this work is only to show that a precise angle 


measuring device is feasible. 


B. LENS TO DETECTOR RAY TRACE 
- This section describes the mathematics of the ray trace 
from the lens to the mirror and then to the detector. A 
listing of program DETPLOT which accomplishes this task is 
included as Appendix C. 
The ray trace consists of the following parts: 
a. Yay direction exiting lens 
b. point of intersection on mirror 
Cc. Yay direction exiting mirror 
d. detector ray count 
1. Ray Direction Exiting Lens 
The ray direction exiting the lens is derived from 
the direction cosines of the ray as calculated py prograin 
GISL. This data is used as an input to program DETPLOT 
along with the coordinates of the point at which the ray 


left the lens. 


2. Point of Intersection on Mirror 
The point of intersection of the ray and mirror may ; 


be described by the following expressions: 


Xm = DemK + XO (5 el) 
Ym = D-mb + Yo (5.2) 
z, = DM + Z5 (5.3) 
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Figure 5.2 Detector Signal. 
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Since the surface of the mirror is tilted along the 

z-axis it can be described by the equations: 
Ym = tan3y(Xm - MP) (5.4) 
Ym cos2 Im + Zone < R (5.5) 

Where equation 5.4 describes a tilted plane and 
equation 5.5 limits the plane to the mirror diameter. 

Neglecting equation 5.5 for a moment, equations 
5.1, 5.2 and 5.4 may be combined to solve for De*,. The 
result is: 
tan Bn Xo ~- MP) - Yo 

gm L- K tan Oe 

Now that Dy, is known, equations 5.1, 5.2 and 5.3 
may be used to find the coordinates of the mirror 
intersection point. These are then checked with equation 
5.5 to make sure they are not outside the mirror diameter. 

3. Ray Direction Exiting Mirror 

The solution of the ray direction leaving the 
mirror is best understood vectorially. See Figure 5.3. 

Reflection off a mirror is governed by two factors. 
Firstly, the angles of incidence and reflection must equal. 


Using vectors this is expressed as: 
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Secondly, the ray is reflected in the same plane as 
the incident ray. This may be expressed vectorially by 
using the cross product to generate the vector Np as 


shown in Figure 5.3. Thus: 


In component form: 
> > > 
ex(—Nyt zi) + ey(-Nyrzi) + ezg(Nyryi-Nyryi) = 


éx (- nyt zy) + ey (“Nyt zr) + ez (Nyrxr ~ Aylyr) 
Three equations result: 
Payee. = Ay Ese 
Shiga See ee 
Mylyi — Aylyi = Ayl yr - Nylyr (5.7) 
The first two results are not unexpected as the 
mirror should not effect the z-component of the ray. 
Equations 5.6 and 5.7 may now be combined to find 
the x- and y-components of the ray. When this is done the 


two following expressions result: 


2 2 
n - on 2n on 
es ake Se a ae 
xv (n2 ¥ n- ) x1 (n2 i n2 ; yl 
x y 4 
(n aoe n 25 2n on 
xr = Yo - i Eig « 
yr (n@ Fe n2 ) yl (n2 ; n2 ) x1 
4 x 4 
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Figure 5.3 Reflection of Ray on Mirror. 
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Since the vector n is entirely in the x-y plane, 
its components may be related to mirror angle as follows: 
Ny = -sin q 
Ny = COS mp 
when the above equations, along with well known 
trigonometric identities are used the following expressions 


will be found. 


Yer cos(2 m)ryji + Sin(2 p)ryi 
yr = -cos(2 m)tyi + Sin(2 m)rxi 
Thus these equations represent the change in direction for 
the exiting ray and form two of its direction cosines, the 
third being frgy-. 
4. Ray Count on Detector 
The coordinates of the ray at the detector are 
described by the expressions: 
Xq = Dng K' + Xm 
Yd = Dna Lt + Ym 
Z2q = Dng M' + 2m 
Since xq is used to define the position of the 
detector, it is known. Thus the value Dpg may be solved 
and the corresponding y and zcomponents inay be found. 
The ray is then checked to make sure it is within 


the detector diameter. If it is the y-coordinate is 


checked to find out if it is negative or positive, and thus 


whether the ray has struck the upper or lower half of the 


detector. The ray count is then incremented accordingly. 
Prograin DETPLOT traces all rays for each increment 
of mirror angle and thus generates plots such as that shown 


in Figure 5.2. 


C. MIRROR AND DETECTOR SYSTEM RESULTS 

Appendix D shows the results of the detector signals as 
a function of mirror angle at angles of incidence from 0 to 
40 degrees. For these the mirror pivot point was placed at 
1.4 units along the x-axis from the inner surface cone 
Since this was the smallest distance that would still allow 
full movement of the scanning mirror. The focal length was 
adjusted to 2.8 units for the lens design calculations so 
that the spot would be imaged as small as possible on the 
detector. Figure 3.4 summarizes the results shown in 
Appendix D. As can be seen, the smaller angles of 
incidence (those in the range 0 to 16 deyrees) result ina 
linear relation between angle of incidence and inirror tilt 
angle, corresponding to an on-axis spot. The laryer 
values of the angle of incidence resulted in a inore random 
relation between these two quantities. 

Appendix D also shows that the error in determining 

mirror angle increases as the angle of incidence increases. 
This is due to the fact that each successive plot is less 


distinct and thus a precise measure of the angle becomes 
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more difficult. Although this effect shows a deterioration 
in the lens system performance, it is not harmful since, as 
the missile corrects itself, the measurement of angle will 
pecome more accurate. 

For these reasons this simple system is considered to 
be a viable means of target detection and shows thnat the 


GISL can be a successful missile seeker lens. 
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Figure 5.4 Mirror Tilt for On-Axis Spot and Angle of 


Incidence, 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Two objectives were accomplished in the course of this 
Study. Firstly, the Gradient Index Seeker Lens design, as 
initiated by Carr [Ref. 4], was completed in that all 
orientations of the center of symmetry with respect to the 
lens were made possible to describe mathematically. The 
results obtained for the spot size on the image plane 
showed trends that indicated the parameters that would 
produce a 'best' lens, but precise definition of these 
values was not possible, 

Secondly, it was shown that a scanning mirror and 
detector system applied to the lens could be described 
using ray-trace techniques and also that such a system 
could be used to provide precise angle measurements of 
off-axis targets at angles less than 16 degrees. For 
larger values of angle, error in definition became 
increasingly greater. 

Specifically, then, the following conclusions are made: 

a) The Gradient Index Seeker Lens can be described 
using ray-trace techniques and can be shown to provide 
adequate imaging of off-axis objects. 

b) A scanning mirror and detector system can be used 


to provide precise angle measurements. 
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B. RECOMMENDATIONS 
The following recommendations for further study are 
made: 


a) An optimization program should be applied to 


program GISL to study the variations of all important 


parameters. 

b) A detector system capable of detection in all 
directions should be designed. 

c) The manner in which signal processing of the 
detector signals can be carried out should be studied. 

dad) A study of Gradient Index material and the methods 


by which a lens could be manufactured should be undertaken. 
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